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The phase diagram of Nf = 3 QCD for smaU baryon densities 
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We demonstrate how to locate the critical endpoint of the QCD phase transition by means of simulations at 
imaginary /x. For the three flavor theory, we present numerical results for the pseudo- critical line as a function of 
chemical potential and bare quark mass, as well as the bare quark mass dependence of the endpoint. 
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1. INTRODUCTION 

As of last year's lattice conference, there were 
three numerical methods ^-[3 capable of simu- 
lating the small fi/T regime of QCD {fJ-B = 3/x) 
with mutually agreeing results on the pseudo- 
critical line Tj, (/is) 0]. However, only one of them 
PP has studied 2-1-1 flavor QCD and obtained a 
prediction for the location of the critical endpoint 
of the deconfinement transition. Here we sum- 
marize new results for three degenerate flavors of 
standard staggered fermions, obtained by simula- 
tions at imaginary chemical potential followed by 
analytic continuation [S] , as previously applied to 
the two flavor case |3] and other observables 
For a discussion of the QCD symmetries at imag- 
inary fi and the general strategy of obtaining the 
critical line we refer to Here we pay special 
attention to the location of the critical endpoint 
as a function of quark mass m. 

In the three-dimensional parameter space 
(T, /i, m) the pseudo-critical temperature repre- 
sents a surface Tc{m,fi). On this surface, the 
line of critical endpoints T*{^) = Tc{mc{lJ,) , fj) 
separates regions of first order transitions from 
crossover behavior. For three standard staggered 
flavors and at zero density, we know the "chiral 
critical point" rndp. = 0) 7 , above which the 
first order deconfinement transition turns into a 
crossover. It belongs to the universality class of 
the 3d Ising model. For this quark mass, the tran- 
sition in the (/i, T)-plane reaches all the way to 
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Figure 1. Critical lines in the (T, /i^)-plane 
for different quark masses, dotted lines denote 
crossover. The bold curve represents T*{fi^). 



the temperature axis, cf. Fig.^ For larger masses 
the critical point moves to positive values of /i, 
while for smaller quark masses it is at imaginary 
values of /i. Our goal is to map out the func- 
tions mc{fJ.),T*{fi) for imaginary fi and analyti- 
cally continue their Taylor expansions to real /i. 

Our simulations were performed on 8^, 10^ and 
12^ X 4 lattices. For each pair (m, /i) we de- 
termined the critical coupling /3c(TO,/i) by in- 
terpolating between three /3-values by means of 
Ferrenberg-Swendsen reweighting For each 
(m, jj,) we accumulated about lOOfc unit length 
trajectories generated by the R-algorithm with 
step size 6t = 0.02. 
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Figure 2. Combined pseudo-critical results for 
L = 8 — 12 and various quark masses. Data for 
different am are shifted to amc(O) according to 
the best fit, which is also shown. 



2. THE PSEUDO-CRITICAL SURFACE 

Data points for all three volumes are shown in 
Fig. 121 They are fitted by truncated Taylor series 



/3c(a/z, am) 



^ CM {a^if^ {am ~ amc{Q))^ . (1) 



Our data are accurate enough to be sensitive to 
next-to-leading terms. A x^-analysis comparing 
all possible fits shows clear preference for a /i^- 
term, while m^- and mixed m/^^-terms appear 
to be smaller than our present errors. On the 
other hand, fits on different volumes are consis- 
tent with each other, so we may fit all volumes to- 
gether. The best fit shown in the plot corresponds 
to Pciatii,am) = 5.1453(2) + 1.780(16)(am - 
0.0323) + 0.705(10)(aAi/)2 + 1.46(15)(aAi/)^. Us- 
ing the two-loop beta-function, this result is easily 
converted to continuum units, 
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-0.602(9) 



1 + 1.937(17) 



me(0) 
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This curve is displayed in Fig.|31together with ear- 
lier results for two |3] and four flavors. Note, 
however, that the present result is the only one 
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Figure 3. Pseudo-critical lines for different 



including a /i^-term. Our three flavor result also 
coincides within its small errors with that for 
Nf = 2-1-1 presented in pQ. This is not sur- 
prising: Tc(/i)/Tc(0) should depend little on the 
quark masses, as long as these are small compared 
to ttT, as is the case in T. 

3. THE LINE OF CRITICAL POINTS 

In order to locate the critical point we used the 
Binder cumulant of the chiral condensate, 



54 = 



(3) 



measured across the pseudo-critical surface. In 
the infinite volume limit, this quantity assumes 
a universal value at a critical point. In the 
universality class of the 3d Ising model it is 
Biimdn), n) = 1.604, while the value is smaller 
for first order transitions and larger for crossover. 
Like for the critical coupling, we fit our data to a 
Taylor expansion 



B4 = 1.604 -I- B (am - amdO) + A{a^iif 



(4) 



Written in this form, A yields directly the desired 
result A = d{amc) / d{aixY . We may now fit 15 
(m, /.t)-pairs measured on our 8^ lattice by a sin- 
gle three-parameter fit, and find A — 0.044(19). 
Data and fit result are shown in Fig. 01 

Unfortunately, the corresponding data for 
10'^, 12'^ are not accurate enough to constrain the 
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parameter A, which is consistent with zero on 
those volumes. However, good fits for B arc ob- 
tained, and those can be used to check consis- 
tency with universahty. The latter predicts that, 
for (3 = Pci'nT'i lAi the correlation length scales as 
^ oc |to - TOc(M)r'', and hence B4 {{L/S,)^/'') = 
1^4 {^L^/'^(am — arricdj.)). We have fitted the vol- 
ume behavior by i? ~ L^^^'^, and obtain v — 
0.62(3) in accord with the 3d Ising value = 0.63. 
Having confirmed universal finite volume scaling 
for our data, we may plot and fit all volumes to- 
gether as in Fig. |S1 The figure nicely demon- 
strates the consistency of our results with uni- 
versal finite size scaling. The fit yields the value 
mc(0) = 0.0323(3), in perfect agreement with ear- 
lier studies done at zero density [7j. However, as 
we mentioned before, the data on the larger vol- 
umes do not constrain the parameter A yet, so 
that we quote our 8'^ number as the final result. 
Converted to continuum units, this yields for the 
/i-dependence of the critical bare quark mass 



rricip- = 0) 



= 1 + 0.84(36) 



ttTJ 



(5) 



This coefficient is inconsistent with a preliminary 
one obtained with improved actions |10| . Action- 
dependent multiplicative renormalization drops 
out of the ratio of quark masses, and any dif- 
ference should be at most 0{a^). Moreover, the 
result violates an analytic upper bound of 9 



(see m). We stress the difficulty of measuring 
the Binder cumulant accurately, without under- 
estimating the error: the tunneling frequency is 
low, and very long MC runs are necessary. 
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